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ABSTRACT 

We present accurate trigonometric parallaxes for 20 new members of the 25 pc white dwarf 
sample as part of the DENSE project (Discovery and Evalution of Nearby Stellar Embers, 
http://www.DenseProject.com). Previously, there were a total of 112 white dwarf systems with 
trigonometric parallaxes placing them within 25 pc and of these, 99 have trigonometric parallaxes 
known to better than 10%. Thus, the 20 new members presented in this work represent a 20% in- 
crease in the number of white dwarfs accurately known to be within 25 pc. In addition, we present 
updated parallaxes for seven known white dwarfs within 10 pc that have been observed as part of 
the ASPENS initiative (Astrometric Search for Planets Encircling Nearby Stars) to monitor nearby 
southern red and white dwarfs for astrometric perturbations from unseen companions. Including a 
few white dwarf companions and white dwarfs beyond 25 pc, we present a total of 33 trigonometric 
parallaxes. We perform atmospheric modeling for white dwarfs to determine physical parameters (i.e., 
Tcff , log g, mass, and white dwarf age). Finally, a new ZZ Ceti pulsating white dwarf was identified 
and revised constraints are placed on two mixed H/He atmosphere cool white dwarfs that display 
continuum absorption in the near-infrared. 

Subject headings: astrometry — Galaxy: evolution — solar neighborhood — stars: distances — white 
dwarfs 
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1. INTRODUCTION 

Knowledge of the low luminosity members of the 
Galaxy — the red dwarfs, subdwarfs, brown dwarfs, and 
white dwarfs (WDs) — comes largely from the nearest 
representatives of each class because they are the eas- 
iest to study. Magnitude-limited surveys of these ob- 
jects allow us to develop crucial nearby samples that re- 
veal ground truths about stellar populations. The most 
straightforward technique to confirm proximity, to amass 
population statistics, and to better constrain physical pa- 
rameters is precision astrometry, specifically trigonomet- 

Electronic address: subasavage@chara.gsu.edu 
Electronic address: jao@chara.gsu.edu 
Electronic address: thenry@chara.gsu.edu 
Electronic address: bergeron@astro.umontreal.ca 
Electronic address: dufourpa@as.arizona.edu 
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ric parallax determinations such as those presented in 
this paper. 

The subjects of this paper, WDs, are perhaps the most 
reliable chronometers of Galactic history available to as- 
tronomers. They also are valuable tracers of Galactic 
populations of different ages, e.g., the thin disk, thick 
disk, and halo. In addition, given that ~95% of all stars 
will end their lives as WDs, the amount of mass contained 
within WDs as a population may already be substantial 
and will continue to increase, possibly comprising a sig- 
nificant fraction of the current (and future) missing mass 
in the Galaxy. 

Here we report results from our effort to enrich the 
sample of WDs within 25 pc by discovering new WDs and 
measuring accurate trigonometric parallaxes for those 
not yet measured. In this paper, we also model nearby 
WDs to determine accurate physical parameters, as well 
as evaluate population statistics for WDs in the So- 
lar neighborhood. Spectral signatures and photomet- 
ric spectral energy distributions (SEDs) are reproduced 
remarkably well by model atmospheres for most WDs. 
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TABLE 1 

Contributions to the 25 pc WD Sample. 



Parallax Program 



All TT TTerr < 10% Refs. 



Yale Parallax Catalog 


104 


91 


1 


Hipparcos 


5 


5 


2,3,4,5 


Torino Observatory Parallax Program 


2 


2 


6 


Ducourant and Collaborators 


1 


1 


7 


CTIOPI 


20 


20 


8 


Total 


132 


119 





References. — (1) van Altena et al. 1995, (2) van Leeuwen 2007, (3) Gould & Chanamc 2004, 

(4) Mugrauer & Neuhauser 2005, (5) Chauvin et al. 2006, (6) Smart et al. 2003, (7) Ducourant et al. 2007, 

(8) this work. 



However, a few exceptional WDs, particularly the coolest 
members, remain problematic. Here we present accurate 
parallaxes for two very cool and low luminosity WDs in 
the solar neighborhood to provide empirical test cases for 
advancements in atmospheric modeling. 

2. CURRENT NEARBY WHITE DWARF CENSUS 

Previous trigonometric parallax efforts have cata- 
logued a total of 112^ WD systems within 25 pc, using no 
uncertainty constraints. A WD system is defined as any 
single or multiple stellar system containing at least one 
WD. The majority of these systems arc included in the 
Yale Parallax Catalog (van Altena et al. 1995) with re- 
cent additions from the Hipparcos space astrometry mis- 
sion (van Leeuwen 2007) as well as ground-based parallax 
programs and companion searches (see Table 1). A com- 
prehensive table containing the current 25 pc WD sam- 
ple has been compiled as part of the DENSE (Discovery 
and Evaluation of Nearby Stellar Embers) project and 
can be found at http://www.DenseProject.com. Surpris- 
ingly, only three WD parallaxes within 25 pc have been 
measured and published since the Hipparcos catalog was 
first released more than a decade ago (ESA 1997)'^. 

To ensure a reliable nearby WD sample, we have 
adopted the quality limit for inclusion into the sample 
that the trigonometric parallax error cannot be larger 
than 10% of the parallax. At 25 pc, this limit amounts to 
an error of 4.0 mas. Given the ~2 mas or better precision 
of ground-based parallaxes, this limit is entirely reason- 
able for the 25 pc sample. Applying this constraint, 13 
systems are eliminated from the 25 pc WD sample, thus 
setting the total number of WD systems with reliable 
parallaxes within 25 pc to 99. 

In order to gauge the degree of incompleteness for the 
25 pc WD sample, we make two basic assumptions, (1) 
that the 10 pc WD sample is complete, and (2) that 
the local density of WDs out to 25 pc is constant. As 
is evident in Figure 1, just over one-third of the WD 
systems expected within 25 pc have trigonometric par- 

^ There are two systems within 25 pc that have radial veloc- 
ity variations which give rise to secondary masses consistent with 
WDs. These systems are G203-047AB (Delfosse et al. 1999) and 
Rcgulus (Gies et al. 2008). We have thus far omitted these objects 
in the statistics pending confirmation. 

^ see note added in manuscript at the conclusion of this article. 



allaxes placing them within that volume. A staggering 
~63% are missing and this is merely a lower limit. If 

additional WDs are found within 10 pc (such as the two 
presented here), the constant density curve shifts verti- 
cally upwards and increases the number of WDs expected 
within 25 pc. The small number of known WDs within 10 
pc (17 systems) presents a fairly large uncertainty when 
extrapolating out to 25 pc solely because of small num- 
ber statistics. However, the fact remains that the sam- 
ple is significantly incomplete. While there are a number 
of known WDs likely within 25 pc that do not yet have 
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Fig. 1. — Cumulative distribution plot for WD systems assuming 
that all WDs out to 10 pc are known (17 systems) and that the local 
WD density is constant out to 25 pc. The solid curve represents the 
number of WD systems expected within a given distance assuming 
constant density. The filled circles are WD systems with accurate 
trigonometric parallaxes within 25 pc. The vertical dashed line 
represents the 25 pc limit. The number of WDs expected (266) vs. 
known (99) within 25 pc is listed to the right of the 25 pc limit. 
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Fig. 2. — Sky distribution plot for the 25 pc WD sample of 99 
systems that meets the 10% parallax error or better criterion {filled 
triangles) as well as the 20 new members of the 25 pc WD sample 
presented here {filled stars) . The overplotted thin lines are Galactic 
coordinate gridlines in increments of 10° with the north and south 
Galactic poles labeled as "NGP" and "SGP", respectively The 
thick line is the Galactic plane with the Galactic center and the 
Galactic anticenter labeled as "GO" and "GAG" , respectively. The 
encircled dot is the direction of the apex of solar motion and the 
encircled cross is the direction of the antapex of solar motion. 

trigonometric parallaxes (e.g., Holbcrg ct al. 2008), there 
remains the need for a sizable sample of nearby, as yet 
undiscovered, WDs to close the incompleteness gap. 

The sky distribution of this sample of 99 WD systems 
is fairly homogeneous (see Figure 2). However, if the 
sky is divided into four equal-area segments by declina- 
tion, it is clear that the southern hemisphere is signifi- 
cantly undersampled in terms of nearby WD systems - 
as shown in Table 2 there are 63 systems in the north 
vs. 36 systems in the south prior to the inclusion of 
the results presented here. In fact, the southernmost 
region is ^50% complete when compared to the north- 
ernmost region and, given that the 25 pc WD sample 
in the northernmost region is likely not complete, this 
is an upper limit. Thus, the southern hemisphere pro- 
vides rich hunting grounds to identify new nearby WDs 
for which our southern hemisphere parallax program, the 
Cerro Tololo Inter- American Observatory Parallax Inves- 
tigation (CTIOPI), is ideally suited. 

In general, there is an observational bias in that higher 
proper motion objects are scrutinized first because they 
are more likely to have larger parallaxes. Thus, there 
may exist a substantial number of WDs with low proper 
motions that remains to be discovered, which may com- 
prise a significant fraction of missing nearby WD sys- 
tems. In order to evaluate empirically whether this 
proposition is plausible, we compared this sample to a 
subset of the Hipparcos catalog. Our intention was to 
evaluate the proper motion distribution of a complete 
sample within 25 pc. Given that this catalog is com- 
plete to F ~ 7.3 - 9.0 (depending on Galactic latitude 
and spectral type. Ferryman et al. 1997), we targeted 
all stars within 25 pc later than AO (assuming My = 
0.6, Binney & Merrifield 1998) down to an apparcnit V 
magnitude of 7.5 {My = 5.5, hence late G type dwarfs). 
A constraint was used to remove evolved stars from the 
sample by defining a line ~1 mag above the locus of 
the main sequence in the H-R diagram and excluding 
those objects with brighter absolute magnitudes. The 



TABLE 2 
25 PC WD Sky Distribution. 



Declination 


# of 


# of New 


Range 


Systems 


Systems 


+90° to +30° 


30 





+30° to +00° 


33 





Total (North) 




63 


-00° to -30° 


21 


7 


-30° to -90° 


15 


13 


Total (South) 




56 



resulting sift included a total of 325 stars, of which 51 
(~16%) had proper motions less than 0'.'15 yr~^, which 
is the detection limit of Lcpine & Shara (2005) and well 
below the detection limits of other proper motion sur- 
veys (e.g., Luyten 1979; Scholz ct al. 2000; Fokorny et 
al. 2003; Hambly et al. 2004)". Thus, it seems the 25 pc 
WD sample is deficient in members with low proper mo- 
tions because only 5 of the 99 systems known before this 
paper have ^ < 0'.'15 yr~^. Age effects, leading to larger 
velocity dispersions for older populations (e.g., WDs), 
are not taken into account that could potentially alter 
these statistics, but likely not by the factor of roughly 
three necessary to resolve the discrepancy. 

In addition to measuring the first trigonometric par- 
allaxes for several nearby WDs, we are conducting an 
effort known as the Astrometric Search for Planets En- 
circhng Nearby Stars (ASPENS, Koerner et al. 2003). 
The ASPENS effort monitors most red dwarfs within 10 
pc and white dwarfs within 15 pc in the southern hemi- 
sphere to search for astrometric signals indicative of un- 
seen companions. The data are acquired in exactly the 
same manner as the parallax data except over longer time 
spans and perhaps with increased cadence. In addition 
to probing for astrometric wobbles, these data allow us 
to redetermine trigonometric parallax;es for nearby WDs 
that are often significantly more accurate than were pre- 
viously available. Ironically, while the southern hemi- 
sphere is undersampled with respect to the 25 pc WD 
sample, the majority of 10 pc WD systems are in the 
south (11 of 17). Thus, we present updated parallaxes 
for seven^ previously known WDs within 10 pc. 

3. OBSERVATIONS AND REDUCTIONS 

Observations have been collected during the ongoing 
CTIOPI program that began in August 1999. CTIOPI 
was conducted first as an NOAO surveys program 
through January 2003 using both the 0.9m and the 1.5m 
(Costa et al. 2005, 2006) telescopes and has since op- 
erated as part of the SMARTS (Smah and Moderate 
Aperture Research Telescope System) Consortium using 

* Deacon & Hambly (2007) conducted an infrared proper motion 
survey down to O'.'l yr~^; however, they employed color constraints 
in search of late-type dwarfs that would eliminate WDs from their 

saiiiplc. 

The remainder of the WD systems in the southern hemisphere, 
except for one - WD 0435-088 (to be published in a future parallax 
publication), are companions to bright stars (e.g., Sirius B) such 
that parallax measurements are not being done on the GTIO 0.9 
m. 
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TABLE 3 
Photometric Results. 



WD Alternate No. of cr No. of No. of 

Name Name Vj .^^KC ^KC Nights tt Filter (mag) Night.s Frames ,7 ctj H (Tjj K S ^o^es 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) 



New 25 pc White Dwarfs 



0038-226 


LHS 1126 


14.50 


14.08 


13.66 


3 


B 





.006 


21 


90 


13.34 





.03 


13,48 





,03 


13 


74 


0-04 




0121-429 


LHS 1243 


14.83 


14.52 


14.19 


4 


R 





.005 


11 


60 


13,86 





,02 


13,63 





,04 


13 


53 


0,04 




0141-675 


LHS 145 


13.82 


13.52 


13,23 


3 


V 





.007 


29 


166 


12,87 





,02 


12,66 





,03 


12 


r,s 


0,03 




0419-487 


GJ 2034 


14.37 


13.76 


12.46 


3 


R 





.010 


12 


64 


10.72 





,02 


10,15 





,02 


9 


85 


0.03 




0628-020 


LP 600 — 42 


15.32 


15.06 


14.75 


3 


I 


[0 


.040] 


15 


70 




















b,c 




LP 600-43 


15.50 


14.13 


12.41 


3 


I 





.007 


15 


74 


10.73 





,03 


10.14 





,03 


9 


,86 


0.02 


o,d 


0751-252 


SCR 0753-2524 


16.27 


15.78 


15.31 


4 


R 


[0 


.029] 


13 


45 


14.75 





,03 


14.47 





,03 


14 


30 


0.09 


b 


0806-661 


L97-3 


13.73 


13.66 


13.61 


4 


R 





.007 


13 


65 


13.70 





,02 


13.74 





,03 


13 


78 


0.04 




0821-669 


SCR 0821-6703 


15.34 


14.82 


14.32 


3 


R 





.007 


17 


86 


13.79 





,03 


13.57 





,03 


13 


,34 


0.04 




1009-184 


WT 1759 


15.44 


15.18 


14.91 


3 


I 





.007 


18 


77 


14.68 





,04 


14.52 





,06 


14 


,31 


0.07 




1036-204 


LHS 2293 


16.24 


15.54 


15.34 


3 


R 





.006 


14 


52 


14.63 





,03 


14.35 





,04 


14 


,04 


0.07 




1202-232 


LP 852 — 7 


12.80 


12.66 


12.52 


3 


R 





.007 


16 


75 


12.40 





,02 


12.30 





,03 


12 


,34 


0.03 




1223-659 


GJ 2092 


14.02 


13.82 


13.62 


3 


V 


[0 


.028] 


14 


61 


13.33 





,04 


13.26 





,06 


13 


,30 


0.06 


b 


1315-781 


L40-116 


16.15 


15.74 


15.36 


3 


R 





.009 


15 


63 


14.89 





,04 


14.67 





,08 


14 


58 


0.12 




1436-781 


LTT 5814 


16.10 


15.81 


15.48 


3 


R 





,006 


18 


63 


15.04 





,04 


14.88 





,08 


14 


76 


0.14 




2008-600 


SCR 2012-5956 


15.84 


15.40 


14.99 


4 


V 





.008 


21 


84 


14.93 





,05 


15.23 





,11 


15 


,41 


Null 




2008-799 


SCR 2016-7945 


16.35 


15.96 


15.57 


4 


R 


0, 


.008 


12 


55 


15.11 





,04 


15.03 





,08 


14 


,64 


0.09 




2040-392 


L495 — 82 


13.75 


13.76 


13.69 


3 


R 





.019 


16 


67 


13.78 





,02 


13.82 





,03 


13 


,81 


0.05 


a,e 


2138-332 


L570-26 


14.48 


14.31 


14.16 


4 


V 





,007 


16 


67 


14.17 





,03 


14.08 





,04 


13 


95 


0.06 




2336-079 


GD 1212 


13.28 


13.27 


13.24 


4 


R 





,009 


19 


74 


13.34 





,03 


13.34 





,02 


13 


35 


0.03 




2351-335 


LHS 4040 


14.52 


14.38 


14.19 


3 


I 


[0 


,043] 


13 


62 


13.99 





,11 


13.86 





,25 


13 


,73 


0.11 


b,c 




LHS 4039 


13.46 


12.33 


10.86 


3 


I 


0, 


,008 


13 


62 


9.48 





,02 


8.91 





,02 


8 


,61 


0.02 


o,f 



Beyond 25 pc White Dwarfs 



0928-713 L64-40 15.11 14.97 

1647-327 LHS 3245 16.20 15.85 

2007-219 GJ 781.3 14.40 14.33 



14.83 3 R 0.006 16 

15.49 3 R 0.009 12 

14.25 3 V 0.007 32 



66 14.77 0.03 14.69 0.06 14.68 0.09 
46 15.15 0.05 14.82 0.08 14,76 0.11 
146 14.19 0.02 14.20 0.04 14.26 0.08 



Known 10 pc White Dwarfs (ASPENS Targets) 



0552-041 LHS 32 14,47 13,99 13,51 3 R 0,006 23 156 13,05 0,03 12,86 0,03 12,78 0.03 

0738-172 LHS 235 13.06 12.89 12.72 4 / 0.009 15 92 12.65 0.02 12.61 0.03 12.58 0.04 

LHS 234 16.69 14.69 12.41 4 / 0.006 15 92 10.16 0.02 9.63 0.02 9.29 0.02 S 

0752-676 LHS 34 13.96 13.58 13.20 3 R 0.006 12 70 12.73 0.02 12.48 0.03 12.36 0.02 

0839-327 LHS 253 11.86 11.77 11.65 3 V 0.007 16 94 11.58 0.03 11.54 0.03 11.55 0.03 

1142-645 LHS 43 11.50 11.34 11.20 3 V 0.007 28 173 11.18 0.01 11.13 0.04 11.10 0.03 

2251-070 LHS 69 15.70 15.11 14.56 3 R 0.007 16 83 14.01 0.03 13.69 0.04 13.55 0.05 

2359-434 LHS 1005 12.97 12.82 12.66 3 R 0.007 12 87 12.60 0.03 12.43 0.02 12.45 0.02 



Likely variable at the --1-2% level (tjee § 5,2), 
^Variability analysis contaminated by nearby source, hence the brackets in column (8) indicating erroneous variability. 

Optical photometry was extracted using PSF fitting rather than by an aperture. 
"^Common proper motion companion to WD 0628 — 020. 
"New ZZ Ceti pulsating WD. 
f 

Common proper motion companion to WD 2351 — 335, 
^Common proper motion companion to WD 0738—172. 



the 0.9m telescope (Jao et al. 2005; Henry et al. 2006). 
The data, results, and procedures presented here corre- 
spond to the 0.9m telescope program. On average, ~80 
nights per year have been allocated to CTIOPI observa- 
tions. The standard CCD setup for CTIOPI observations 
(both photometric and astrometric) utilizes only the cen- 
tral quarter of the 2048 x 2046 Tektronix CCD camera 
with 0'.'401 pixel^^, yielding a 6.'8 square field of view. 
The Tek 2 Vj-RKC-^KC^ (hereafter without the subscripts) 
filter set was used to carry out the observations. 

3.1. Photometry 

Photometric observations have been collected since the 
inception of CTIOPI during scheduled observing runs 
when sky conditions were photometric. Standard stars 

from Graham (1982) and Landoh (1992, 2007) were 
taken nightly through a range of airmasses to calibrate 

^ The central wavelengths for Vj, i?KCi ^KC ^I'e 5475, 6425, 
and 8075 A, respeetively. The Tek 2 Vj filter cracked in 2005 
March and was replaced by the very similar Tek 1 Vj filter. See § 
3.2.3 for a discussion on the impact this switch has on the data. 



fluxes to the Johnson-Kron- Cousins system and to cal- 
culate extinction corrections. Bias subtraction and fiat- 
fielding (using calibration frames taken nightly) were per- 
formed using standard IRAF packages. An aperture of 
14" diameter was used when possible (consistent with 
Landolt 1992) to determine stellar flux. Cosmic rays 
within this aperture were removed before flux extrac- 
tion. Aperture corrections were applied when neighbor- 
ing sources fell within the adopted aperture. In these 
cases, the largest aperture that did not include flux from 
the contaminating source was used and ranged from 4" 
to 12" in diameter. Total uncertainties (including inter- 
nal night-to-night variations as well as external flts to 
the standard stars) in the optical photometry are ± 0.03 
mag in each fllter (Henry et al. 2004). In the cases of 
WD 0628-020 and WD 2351-335, the primaries (red 
dwarfs) were a significant contaminant (especially in /) 
so that aperture photometry alone was not possible. In- 
stead, a PSF fit using interactive data language (IDL) 
package mpfit2dpeak was generated for the primaries and 
then subtracted from the images. Aperture photometry 
was then performed on the subtracted images. The same 
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procedure was performed for extracting photometry for 
the primaries (i.e., the WD secondaries were PSF fit and 
removed). 

Our photometric results are given in Table 3, which 
is divided into three samples based on the trigonometric 
parallaxes presented here: (1) new 25 pc WD members, 
(2) WDs beyond 25 pc, and (3) known 10 pc ASPENS 
targets. Companions to WDs (i.e., LP600-43, LHS 234 
and LHS 4039) for which photometric analyses were per- 
formed are included in the table just below their WD 
companions. 

Multi-epoch optical VRI photometry [columns (3)- 
(5)] as well as near-infrared JHKs photometry [columns 
(11)-(16)] extracted from the 2MASS database are Hsted 
in Table 3. We performed a photometric variability anal- 
ysis of the parallax target (hereafter referred to as the 
"PI star") relative to the reference stars using the par- 
allax data taken in a single filter, as outlined in Honey- 
cutt (1992). Columns (6)-(10) fist the number of differ- 
ent nights the object was observed for photometry, the 
parallax filter, the standard deviation of the PI star's 
magnitude in that filter from the variability analysis, the 
number of nights parallax data were taken, and the total 
number of parallax frames used in the variability anal- 
ysis. In general, objects with standard deviation val- 
ues larger than 0.02 mag are considered variable, those 
with values between 0.01-0.02 mag are likely variable at 
a few percent level, and those less than 0.01 mag are 
"steady". There are four cases where the standard de- 
viations are larger than 0.02 mag (WD 0628-020, WD 
0751-252, WD 1223-659, and WD 2351-335) and in 
all four cases, there are contaminating sources within a 
few arcseconds. Thus, the large standard deviations are 
due to varying degrees of contamination within the pho- 
tometric aperture depending on seeing conditions rather 
than intrinsic variability. Two objects show standard de- 
viations between 0.01-0.02 mag (WD 0419-487 and WD 
2040-392). WD 0419-487 has an unresolved red dwarf 
companion in a short period orbit (see § 5.2) so that mild 
variability might be expected. WD 2040—392 is a new 
pulsating ZZ Ceti WD (see § 5.2). 

3.2. Astrometry 

A complete discussion of parallax data acquisition and 
reduction techniques can be found in Jao et al. (2005). 
Briefly, once an object (or system) is selected to be ob- 
served for a trigonometric parallax determination, the 
object needs to be "set up" . This process consists of se- 
lecting a telescope pointing as well as a filter bandpass 
{VRI) with which all subsequent astrometry observa- 
tions will be taken. The telescope pointing is selected so 
that a fairly homogeneous distribution of reference stars 
encircle the PI star, do not reside near bad columns, and 
are as close as possible to the PI star rather than near 
the edges of the CCD. Also, every efli'ort is made to place 
the PI star as close to the center of the CCD as possi- 
ble. Of course, compromises are made in cases of sparse 
fields. 

Once sufficient data have been collected (criteria that 
define a definitive parallax are outlined in § 3.2.1), each 
frame is inspected and poor quality frames (e.g., bad 
seeing, telescope guiding problems) are discarded. Cen- 
troids for the reference field and PI star are extracted 
using SExtractor (see discussion in § 3.2.2, Bertin & 



Arnouts 1996). The centroids are corrected for differen- 
tial color refraction (DCR) based on the color of each ref- 
erence star and PI star. Typically, this correction shifts 
the stars' positions by no more than a few mas. 

Automatic quality control constraints on ellipticity, 
elongation, and full width at have maximum (FWHM) 
eliminate frames and individual stars (reference and PI) 
of poor quality not noticed during manual inspection. 
A fundamental trail plate is chosen as a reference plate 
and rotated based on a comparison with either the Guide 
Star Catalog 2.2 or 2MASS (used by default except when 
there are too few stars in the PI star field that are cat- 
alogued in the 2MASS database). All star centroids on 
other plates are recalibrated to account for different scal- 
ing in both the X and Y directions, as well as the different 
amounts of translation and rotation. A least-squares re- 
duction via the Gaussfit program (.Jefferys et al. 1988) is 
performed, assuming the reference star grid has STTi = 
and S/ij = 0, where tt and are parallax and proper mo- 
tion, respectively. After verification of a good reference 
field (e.g., no reference stars within 100 pc, no problem- 
atic reference stars because of unresolved multiplicity), 
the final reduction produces a relative trigonometric par- 
allax for the PI star. 

Even the reference stars trace out small parallax el- 
lipses because these stars are not infinitely far away. 
Thus, a correction to absolute parallax must be per- 
formed. This is accomplished in one of (at least) three 
different ways, (1) using a model of the Galaxy for the 
disk and halo, (2) spectroscopic parallaxes for each of the 
reference stars, or (3) photometric parallaxes for each 
of the reference stars. Because accurate photometry is 
needed for DCR correction and is already available, we 
use photometric parallaxes for the reference stars to cor- 
rect to absolute parallax using the CCD distance rela- 
tions of Henry et al. (2004). These relations assume the 
reference stars are single main-sequence dwarfs and do 
not take into account contamination from evolved stars, 
unresolved double stars, or reddening. We identify and 
remove outlying points estimated to be within ~100 pc 
due either to evolved stars or unresolved double stars 
whose distances are underestimated. Given that all of the 
WDs presented here are relatively nearby, PI star red- 
dening is negligible. However, the distant reference stars 
can be reddened and appear to be closer than they truly 
are. This is the case for two WD fields (WD 1223-659 
and WD 1647—327) for which we have adopted an aver- 
age correction to absolute of 1.2 mas (based on all other 
WD fields presented here) with a conservative error of 
0.3 mas (consistent with the largest absolute correction 
errors presented here). 

3.2.1. Definitive Parallax Criteria 

During the course of CTIOPI, hundreds of parallaxes 
(both preliminary and definitive) have been reduced. 
Based on our experience with the "evolution" of the 
parallaxes as more data are added, we have set limits 
that define when a reduction is sufficient to be deemed 
definitive and hence publishable. In order to accurately 
decouple the parallactic and proper motions in the fi- 
nal astrometric solution, observations must span at least 
two years and adequately sample the parallax ellipse (in- 
cluding high parallax factor observations). Also, there 
must be a balance between negative parallax factor ob- 
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servations and positive parallax factor observations and 
typically at least 20 good frames of each are required. At 
least two and ideally three VRI photometry observations 
per PI star are necessary to ensure accurate photometry 
needed to correct for DCR as well as to absolute paral- 
lax. Finally, the parallax error must be less than 3.0 mas. 
This constraint was set early on during the program and 
with the use of the new centroiding algorithm (see next 
section) is typically an easy requirement to meet. 

3.2.2. New SExtractor Centroiding Algorithm 

Beginning with SExtractor version 2.4, the authors im- 
plemented windowed positional parameters. These serve 
to alleviate a number of the inherent inaccuracies of the 
isophotal positional parameters used for centroiding in 
previous versions of Sextractor (see S extractor y2.0 User's 
Guide for a complete discussion). The authors analo- 
gize the positional accuracy of the Gaussian-weighted 
two-dimensional windowed centroids to that offered by 
point-spread-function (PSF) fitting. After extensive test- 
ing on dozens of PI stars from the CTIOPI program, 
we found that windowed centroids were clearly superior 
to the isophotal centroids used in our previous reduc- 
tions, including those published in Jao et al. (2005) and 
Henry et al. (2006) from the 0.9m. On average, paral- 
lax errors were reduced by ~50% with no appreciable 
change in the parallax (when PI star contamination was 
not present) using the same data sets in both reductions. 
In extreme cases, particularly in crowded fields, parallax 
errors were reduced by factors of two or more. This was 
the case with the parallax reduction for WD 0751—252, 
which had a contaminating source within I'.'O. Using the 
isophotal centroids, the preliminary parallax determina- 
tion was 68.6 ± 3.5 mas with clear trends in the residuals 
indicative of contamination. Using the windowed cen- 
troids and the same data set, the parallax determination 
was 55.8 ± 1.0 mas with residuals that hovered near zero 
throughout the data^. As an additional check, this WD 
was discovered to be a common proper motion compan- 
ion to LTT 2976 (Subasavage et al. 2005b), which has 
a Hipparcos parallax of 51.52 ± 1.46 mas (van Leeuwen 
2007). While this value deviates by slightly more than 
2(7 of the value determined for WD 0751—252, it is in far 
better agreement than the result using the isophotal cen- 
troids. We obtain an average parallax error of 1.10 mas 
for the sample presented here using windowed centroids 
(excluding WD 0628-020 and WD 2351-335 because 
the observations were optimized for their red dwarf com- 
panions and as such, the exposure times were shorter and 
the WD components were underexposed). 

3.2.3. Cracked V Filter 

In early 2005, the standard Tek 2 V filter that was 
included in the VRI filter set used for parallax observa- 
tions cracked. It was replaced with another V filter that 
had a very similar transmission profile. Upon reduction 
of parallaxes including data taken in mid-2005, the tar- 
gets that were observed in the V filter appeared to show 
subtle arches in the residuals once parallax and proper 

^ Recent data not used in the centroiding test as well as a correc- 
tion to absolute parallax are included in the definitive parallax of 
56.54 it 0.95 mas presented here, hence the slightly different value 
in Table 5. 
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Fig. 3. — Hertzsprung-Russel diagram with the new WD systems 
presented here overplotted on a sample of known WDs within 25 
pc from Bergeron et al. (2001). The system labeled 0419-487 is 
discussed in § 5.2. 

motion were fit and removed. It became more evident 
as additional data were collected. Now that we have ~3 
years of data after the filter switch, we are able to see its 
effect as a "dip" in the PI star residuals of a few mas for 
data taken in mid-2005. We have performed extensive 
tests and have concluded that inclusion of both V filters' 
data in the final parallax reduction is beneficial for con- 
straining the parallax provided there are sufficient data 
that span at least ~l-2 years after the filter switch. The 
major drawback is that any astrometric perturbations 
from unseen companions will likely be missed (unless the 
photocentric shift is large) because the residuals are con- 
taminated. Given that the majority of WD parallaxes 
presented here are taken in R (19 vs. 7 a,t V and 4 at /) 
the effect is minimal for this sample. 

4. ASTROMETRY RESULTS 

Astrometric results for the WD systems (including 
companions for which astrometric analyses were per- 
formed) are listed in Table 5. Columns (4)-(9) list the 
filter used for parallax observations, the number of sea- 
sons the PI star was observed, total number of frames 
used in the parallax reduction, the time coverage and 
length of the parallax data, and the number of reference 
stars used. The 'c' in column (5) signifies that the obser- 
vations were continuous throughout every season within 
the time coverage. The 's' signifies that observations 
were scattered such that there is at least one season with 
only one night's data (or no data for an entire season). 
Columns (10)-(12) list the relative parallax, correction to 
absolute, and the absolute parallax. The proper motions 
and position angles quoted in columns (13) and (14) are 
those measured with respect to the reference field (i.e., 
relative, not corrected to a non-rotating frame of refer- 
ence). The tangential velocities quoted in column (15) 
are not corrected for solar motion. 

As can be seen in the observational Hertzsprung-Russel 
(H-R) diagram in Figure 3, all but one of the observed 
WDs lie within the realm of known WDs. The only ex- 
ception is WD 0419—487 for which an unresolved red 
dwarf companion significantly contaminates the optical 
and near-IR photometry (see § 5.2). In Figure 4, CTIOPI 
parallaxes for WDs previously known to be within 10 pc 
(ASPENS targets) are compared to their parallaxes from 
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Fig. 4. — Comparison plots of the CTIOPI parallaxes vs. pho- 
tographic plate-derived parallaxes from the Yale Parallax Catalog 
(YPC, van Altena et al. 1995) and in the case of WD 1142-645, 
vs. the recently updated Hipparcos parallax (HIP, van Leeuwen 
2007). 

the Yale Parallax Catalog (YPC, van Altena et al. 1995) 
and, in one case (WD 1142—645), with the Hipparcos 
parallax (van Leeuwen 2007). Agreements are reasonably 
good (within 2a), especially considering that all of the 
values from the YPC were obtained using photographic 
plates. 

Of the new 25 pc WD members from this effort, the 
majority tend to have proper motions on the lower end 
of the distribution (see Figure 5 a). Based on the dis- 
cussion in § 2, this trend is to be expected and likely 
represents the realm where the majority of nearby WDs 
as yet undiscovered will be found. For instance, ~90% 
of the 25 pc WD sample members have proper motions 
greater than WD 1202-232 with /x = 0'.'246 yr'^ yet it 
is now the 25th nearest WD system (10.83 ± 0.11 pc), 
and WD 2336—079 has the lowest proper motion of the 
sample with /i = 0'.'034 yr^^. While none of the new 25 
pc members are within 8 pc, five were found between 8 
and 15 pc (see Figure 56). This includes one object with 
a poorly constrained previous parallax (WD 0038—226, 
van Altena ct al. 1995), two that have been suspected to 
be nearby for many years but whose trigonometric paral- 
laxes were not determined until now (WD 0141—675 and 
WD 1036-204, Gliese & Jahreifi 1991), and two recent 
discoveries (WD 0821-669 and WD 1202-232, Subasav- 
age et al. 2007). Of particular interest is that the two 
recent discoveries arc closer than the 13 pc limit set by 
Holberg et al. (2002) from which they determine the local 
WD density. While the most recent local WD density de- 
termination (Holberg et al. 2008) takes these two objects 
into account, they support the idea that more nearby 
WDs may yet be found. The distribution of tangential 
velocities for the new 25 pc members is unexceptional 
(see Figure 5c). Only one object (WD 2008-600) has a 
tangential velocity greater than 100 km sec~^ (see § 5.2). 

5. ANALYSIS 

5.1. Modeling of Physical Parameters 

Atmospheric modeling procedures of the WDs are 
identical to those presented in Subasavage ct al. (2008) 
and references therein, with the exception that the 
trigonometric parallax constrains the surface gravity (in- 
stead of assuming log g = 8 as was done in that publica- 
tion). Briefly, optical/near-IR magnitudes are converted 
into fluxes using the calibration of Holberg et al. (2006) 




1 2 3 4 5 5 10 15 20 25 100 200 300 

/X (arcsec yr"') Dislance„ (pc) Vjgn (km s"') 

Fig. 5. — Histograms of number of WDs within 25 pc vs. (a) 
proper motion binned by 0'.'5 yr^-'^ (6) distance binned by 2 pc and 
(c) tangential velocity binned by 20 km sec In all three his- 
tograms, the unshaded region represents the 99 previously known 
WDs within 25 pc and the shaded region corresponds to the 20 
new discoveries presented in this work. The dip in panel 6 at ~14 
pc is an artifact of the binning. 

and compared to the SEDs predicted by the model at- 
mosphere calculations. The observed flux, /™, is related 
to the model flux by the equation 

f^ = A7r{R/Df (1) 

where R/D is the ratio of the radius of the star to its 
distance from Earth, 7J™ is the Eddington flux (depen- 
dent on Tgff, logg, and atmospheric composition) prop- 
erly averaged over the corresponding fllter bandpass, and 
TT in this context is the mathematical constant (elsewhere 
throughout this paper, tt refers to the trigonometric par- 
allax angle) . Our fitting technique relics on the nonlinear 
least-squares method of Levenberg-Marquardt (Press et 
al. 1992), which is based on a steepest descent method. 
The value of is taken as the sum over all bandpasses 
of the difference between both sides of Equation (1), 
weighted by the corresponding photometric uncertain- 
ties. Only Teff and [Tr{R/D)'^] are free parameters and the 
uncertainties of both parameters are obtained directly 
from the covariance matrix of the fit. The main atmo- 
spheric constituent (hydrogen or helium) is determined 
by the presence of Ha from spectra published in the lit- 
eratuK! (references listed in Table 4) or by comparing fits 
obtained with both compositions. 

We start with log g = 8.0 and determine T^s and 
[tt{R/ D)'^], which combined with the distance D ob- 
tained from the trigonometric parallax measurement 
yields directly the radius of the star R. The radius 
is then converted into mass using evolutionary mod- 
els similar to those described in Fontaine et al. (2001) 
but with C/0 cores, g(He) = logM^e/M^, = 10~^ and 
(jf(H) = 10~^ (representative of hydrogen-atmosphere 
WDs), and q(Rc) = 10"^ ^nd g(H) = 10"^" (repre- 
sentative of helium-atmosphere WDs).* In general, the 
log g value obtained from the inferred mass and radius 
{g = GM/R^) will be different from our initial guess of 
log g = 8.0, and the fitting procedure is thus repeated 
until an internal consistency in log g is reached. The 
parameter uncertainties are obtained by propagating the 
error of the trigonometric parallax; measurements into 
the fitting procedure. 

Physical parameter determinations for the DQ and DZ 
WDs are identical to the procedures outlined in Dufour et 
al. (2005, 2007). Briefly, the photometric SED provides a 
first estimate of the atmospheric parameters with an as- 
sumed value of metal abundances using solar abundance 
ratios. The optical spectrum is fit to better constrain 
the metal abundances and to improve the atmospheric 

^ see http://www.astro.umontreal.ca/~bergeron/CoolingModels/. 
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TABLE 4 
Physical Parameters. 



WD Adopted tt'^ No. Tg£f Age*^ 

Name (mas) of tt Ref. (K) log g Comp Ref. M/Mq My log L/Lq (Gyr) Notes 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 



New 25 PC White Dwarfs 



003S- 


226 


0121- 


429 


0141- 


675 . . 


0419- 


487 . . 


0628- 


020 . . 


0751- 


252 


0806- 


661 , 


0821- 


669 


1009- 


184 


1036- 


204 


1202- 


232 


1223- 


659 , 


1315- 


781 , 


1436- 


781 , . 


2008- 


799 . . 


2008- 


600 . . 


2040- 


392 . . 


2138- 


332 . . 


2336- 


079 . . 


2351- 


335 . . 



110. 


42±1 


17 


2 


1.2 


5210±130 


7.92±0 


.02 


He(+H) 


6 





.52±0. 


.01 


14 


.72±0 


04 


-3 


.94±0 


.01 


4 


.37±0, 


.18 


c 


54, 


,61±0 


.96 


1 


2 


















13 


.52±0, 


.05 














d 


102. 


80±0 


.85 


1 


2 


6470±130 


7.99±0 


.01 


H 


7 





.5S±0 


.01 


13 


.88±0 


.03 


-3 


.58±0 


.01 


1 


.81±0, 


.04 




49, 


.68±1 


.34 


1 


2 


















12 


.85±0, 


.07 














e 


48 


.84±1. 


.10 


2 


2,4 


6846±383 


8.06±0. 


.04 


H 


8 


0, 


.63±0. 


.02 


13 


.76±0. 


.06 


-3 


.53±0 


.02 


1. 


.78±0, 


.10 




55, 


.05±0 


.80 


2 


2,3 


5160±100 


7.97±0 


.02 


H 


9 





.56±0, 


.01 


14 


.97±0 


.04 


-3 


.97±0 


.02 


4 


.46±0 


.30 




52 


17±1 


.67 


1 


2 


10250± 70 


8.06±0 


05 


Ho( + C) 


10 





62±0 


03 


12 


.32±0 


08 


-2 


.83±0 


.03 





67±0 


04 


t 


93 


89±1 


04 


1 


2 


5150±100 


8.11±0 


02 


H 


11 





65±n 


01 


15 


.2n±o 


04 


-4 


.06±0 


.01 


6 


22±0 


16 




55 


.55±0 


.85 


2 


2,3 


5940±280 


8.05±0 


.02 


He( + H,Ca) 


10 





,60±0 


.01 


14 


.16±0 


.04 


-3 


.78±0 


.02 


2 


.97±0 


.18 




70, 


.00±0 


.66 


1 


2 


















15 


.47±0, 


.04 














g 


92, 


.37±0, 


.91 


1 


2 


8590±170 


7.90±0, 


.02 


H 


12 


0, 


,54±0, 


.01 


12 


.63±0, 


.04 


-3 


.04±0 


.01 


0, 


.78±0, 


.02 




61 


,43±1 


16 


1 


2 


7660±220 


7.81±0 


03 


H 


13 





49±0 


02 


12 


.96±0 


05 


-3 


.19±0 


.02 





93±0 


04 




52 


14±0 


94 


1 


2 


5730±160 


8.18i0 


03 


H 


11 





70±0 


02 


14 


.74±0 


05 


-3 


.91±0 


.02 


4 


15±0 


22 




40, 


,56±0 


.80 


1 


2 


6270±200 


8.08±0 


.03 


H 


11 





,64±0 


.02 


14 


.14±0 


.05 


-3 


.70±0 


.02 


2 


,41±0 


.18 




40, 


,06±1, 


35 


1 


2 


5800±160 


7.97±0. 


06 


H 


11 


0, 


.57±0. 


03 


14 


.36±0. 


08 


-3 


.77±0 


.03 


2 


,34±0, 


25 




60 


.42±0. 


86 


1 


2 


50S0±220 


7.89±0. 


03 


He( + H) 


11 





.51±0. 


02 


14 


.75±0. 


04 


-3 


.97±0, 


.02 


4 


,55±0, 


27 


c 


44 


.18±0. 


.97 


1 


2 


10S30±310 


S.03±0. 


.04 


H 


11 





.62±0. 


.02 


11 


.98±0. 


.06 


-2 


.71±0, 


.02 


0, 


.51±0. 


.02 




64, 


.00±1. 


.41 


1 


2 


7240±260 


S.18±0. 


.03 


He(+Ca) 


11 





.69±0. 


.02 


13 


.51±0. 


.06 


-3 


.51±0 


.02 


1 


.96±0. 


.12 




62, 


.72±1. 


.70 


1 


2 


11000±300 


S.25±0. 


.04 


H 


14 





.76±0. 


.02 


12 


.27±0. 


.07 


-2 


.81±0 


.03 





.69±0. 


.04 




43 


.74±1. 


.43 


2 


2,4 


8070±390 


7.80±0. 


.05 


H 


15 





.49±0. 


.03 


12 


.72±0. 


.08 


-3 


.10±0, 


.03 





.S1±0. 


05 





Beyond 25 pc White Dwarfs 



0928-713 38.44±0.64 1 2 

1647-327 37.03±1.18 1 2 

2007-219 38.22±0.94 1 2 



8880±260 8.25±0.03 H 13 

6120±200 7.92±0.05 H 11 

9520±230 7.97±0.04 H 16 



0.75±0.02 13.03±0.05 -3.19±0.02 1.18±0.04 
0.54±0.03 14.04±0.08 -3.64±0.03 1.87±0.14 
0.58±0.02 12.31±0.06 -2.90±0.03 0.66±0.03 



Known 10 pc White Dwarfs (ASPENS Targets) 



0552-041 155.97±0.78 2 1,2 5180± 70 8.35±0.01 He(+H,Ca) 17 0.80±0.01 15.44±0.03 -4.20±0.01 6.82±0.02 

0738-172 109.43±0.56 3 1,2,4 7600±220 8.03±0.01 He( + H,Ca) 17 0.60±0.01 13.26±0.03 -3.34±0.01 1.41±0.02 

0752-676 126.62±1.32 2 1,2 5700± 90 8.00±0.02 H 6 0.59±0.01 14.47±0.04 -3.S1±0.01 2.65±0.10 

0839-327 113. 59± 1.93 2 1,2 9120±190 7.72±0.03 H 18 0.45±0.01 12.14±0.05 -2.84±0.02 0.55±0.02 

1142-645 216.12±1.09 3 1,2,5 7920±220 8.07±0.01 He(+C) 19 0.62±0.01 13.17±0.03 -3.29±0.01 1.32±0.01 

2251-070 117.38±0.95 2 1,2 4000±200 7.92±0.02 He(+H,Ca) 17 0.52±0.01 16.05±0.03 -4.40±0.01 7.39±0.18 

2359-434 122. 70± 1.13 2 1,2 8530±160 8.39±0.01 H 16 0.85±0.01 13.40±0.04 -3.35±0.01 1.82±0.06 



References. — (1) YPC (van Altena et al. 1995) tt, (2) this work tt, (3) Hipparcos companion tt (van Leeuwen 2007), (4) this work companion tt, 
(5) Hipparcos (van Leeuwen 2007) tt, (6) Bergeron et al. 1997, (7) Hintzen &; Jensen 1979, (8) Silvestri et al. 2001, 

(9) Subasavage et al. 2008, (10) this work, (11) Subasavage et al. 2007, (12) Kilkenny et al. 1997, (13) Wickramasinghe & Bcssell 1977, 
(14) Berger Frlngant 1984, (15) McCook & Sion 1999, (16) Eggen &c Greenstein 1965, (17) Dufour et al. 2007, 
(18) Bergeron et al. 2001, (19) Dufour et al. 2005. 

*The adopted parallaxes are weighted means in cases of multiple parallax determinations for a system. Model parameters were determined using these values. 
The Ref. column (4) identifies the source of each parallax. 
^WD cooling age only, not including main-sequence lifetime. 

'^Atmospheric modeling included using trace hydrogen in a helium atmosphere to best reproduce the SED (sec § 5.2). 
*^Physical parameters are not listed because there is evidence that this object is an unresolved double degenerate (see § 5.2). 

^Object is an unresolved red dwarf-white dwarf binary whose photometry and spectroscopy are strongly contaminated so that no atmospheric modeling was 
possible. 

^Atmospheric modeling included the ultraviolet spectrum from lUE as well as the optical/near-IR photometry (see § 5.2). 
^Atmospheric modeling was not possible because the sourcc(s) of spectral features not yet well understood. 
Effective temperature is the limit of the model grid and additional pressure effects in this regime are not accounted for. 



parameters from the photometric SED. This procedure 
is iterated until a self-consistent photometric and spec- 
troscopic solution is reached. 

Only two objects' spectra are modeled here for the first 
time, WD 0806-661 and WD 1009-184. For the remain- 
ing DQ (those with carbon) and DZ (those with calcium) 
stars, spectral modeling to obtain abundances was per- 
formed and presented in Dufour et al. (2005, 2007) and 
Subasavage et al. (2007, 2008). The atmospheric abun- 
dances will not change with the inclusion of the paral- 
laxes; however, the surface gravities (hence masses) are 
sensitive to changes in distance and have been updated in 
Table 4. For WD 0806—661, the optical spectrum shows 
no carbon absorption, yet it is classified as a DQ based 
on ultraviolet (UV) spectra. Thus, the UV spectrum was 
used for fitting (see § 5.2). For DZ stars, trace amounts 
of hydrogen not directly visible can be present in the at- 
mosphere and affect the spectral profiles of the calcium 



absorption lines (Dufour et al. 2007). In the case of WD 
1009—184, whose spectrum was obtained using the same 
telescope/instrument setup and reduction procedures as 
described in Subasavage et al. (2008), the spectral fit 
including a log (H/He) = —3 better reproduced the cal- 
cium H & K lines than if no hydrogen were present (see 
§5.2). 

In order to best constrain the physical parameters for 
this sample, weighted mean parallaxes and errors are 
calculated for systems with previous parallax determi- 
nations as well as those that have common proper mo- 
tion companions with previous/new parallax determina- 
tions. The parallax values that are used to model the 
physical parameters are listed in Table 4 [column (2)] as 
well as the number of individual parallaxes used in the 
weighted mean and corresponding references [columns 
(3) and (4)]. Colmnns (5) and (6) list the effective tem- 
peratures and surface gravities as well as corresponding 
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errors. Columns (7) and (8) list the composition(s) used 
in the atmospheric modeling (with any secondary con- 
stituents listed in parentheses) and the spectral type ref- 
erence. Columns (9)-(13) list the derived masses, abso- 
lute magnitudes, luminosities, WD ages (not including 
main-sequence lifetimes), and any notes for the systems. 
While the errors listed for mass, luminosity, and age are 
formal errors, they are remarkably well constrained when 
accurate trigonometric parallaxes are available. 

5.2. Comments on Individual Systems 

WD 0038-226: this WD has mild absorption bands 
similar to the carbon Swan bands found in DQ stars but 
the bands are shifted blueward. Initially, it was thought 
that these bands were actually the Swan bands but were 
pressure-shifted because of increased pressures in cool 
He-rich WD atmospheres (Liebcrt & Dahn 1983). An- 
other explanation was that the features were caused by 
the hydrocarbon C2H (Schmidt et al. 1995). The most 
recent possible explanations revisit the idea of pressure- 
shifted Swan bands alone or in conjunction with Swan 
bands produced by highly rotationally excited C2 (Hall 
& Maxwell 2008). Additional theoretical investigations 
are necessary to better understand the properties of C2 in 
the high-pressure, high-temperature helium environment 
that a WD atmosphere would provide. We present a sig- 
nificantly better parallax (the previous parallax had an 
error — 10.3%), that confirms this object is within 10 pc 
and is the nearest known WD with the aforementioned 
spectral anomaly. 

In addition. Bergeron et al. (1994) have shown that 
this object displays collision-induced absorption in the 
infrared and can be attributed to collisions of molecu- 
lar hydrogen with helium. Bergeron et al. (1997) have 
modeled the SED using a mixed H/He composition and 
arrived at a satisfactory fit. The updated physical pa- 
rameters found in Table 4 were derived by an identical 
analysis (including the use of their B VRIJHK photome- 
try) except with the updated trigonometric parallax pre- 
sented here. The helium abundance derived with the 
updated parallax [log (He/H) = 1.31] is less than that 
found by Bergeron et al. (1997) [log (He/H) = 1.86]. 

WD 0121—429: a recently discovered magnetic DA 
that is thought to be an unresolved double degenerate 
with one component being a magnetic DA and the other 
being a featureless DC (Subasavage et al. 2007). The 
physical parameters are not listed in Table 4 because 
any number of component masses and luminosities can 
reproduce the SED fit. 

WD 0141-675: a DA WD that was one of only two 
new 25 pc WD members to be within 10 pc. 

WD 0419—487: a WD + red dwarf pre-cataclysmic 
eclipsing binary also known as RR Caeli with an orbital 
period of 7.3 hours for which Maxted et al. (2007) de- 
rive masses and radii. Contamination from the main- 
sequence component is significant both spectroscopically 
and photometrically so that no atmospheric modeling of 
the WD was possible (it is the outlying point in the H-R 
Diagram in Figure 3). 

WD 0628-020: a known WD that has a red dwarf 
companion (exact spectral type unknown) with a separa- 
tion of 4'.'5 at 317.9°. Parallax data were taken in / to op- 
timize observations for the red dwarf companion. Thus, 
the parallax for the WD is more uncertain because it is 



somewhat poorly exposed in our images. The weighted 
mean parallax listed in Table 4, which represents our 
two measurements, should be taken as the distance to 
the system. 

WD 0738-172: a known nearby WD that has a 
M6.0V companion with a separation of 20'.'6 at position 
angle 261.0°. Because the companion is fainter in both 
V and R as well as redder, the parallax observations 
were taken in / so that an independent parallax to the 
secondary red dwarf (of comparable brightness to the 
primary at /) could also be obtained. The parallax val- 
ues for the pair are in excellent agreement (see Table 5) 
but the proper motion values differ by several sigma of 
the formal errors. By comparing astrometric results of 
several distant wide binaries on CTIOPI (as yet unpub- 
lished), we find that the proper motion and position angle 
values agree to within 1-2(T. Thus, the formal errors are 
likely not significantly understated. We make some ba- 
sic assumptions about component masses and orientation 
(i.e., face-on orbit) of the system to evaluate the plau- 
sibility of orbital motion to account for the discrepancy. 
Indeed, orbital motion during four years of observation 
of this nearby (9.14 ± 0.05 pc) system is likely the cause 
of the discrepancies in /i and position angle. 

WD 0751-252: a recently discovered WD that is a 
common proper motion companion to LTT 2976 (Sub- 
asavage et al. 2005b). The trigonometric parallax deter- 
mined in this work for the WD (56.54 ± 0.95 mas) is in 
marginal agreement with the Hipparcos parallax for the 
primary (51.52 ± 1.46 mas, van Leeuwen 2007). 

WD 0806-661: a DQ WD that shows carbon fea- 
tures only in the UV (in the optical, it appears as a 
featureless DC). Thus, a UV spectrum from the Interna- 
tional Ultraviolet Explorer (lUE) archive was included 
in the analysis. First, a model that includes carbon with 
an abundance below the detectability limit in the optical 

[e.g., log (C/He) 4 to -6] was fit to the VRIJHK 

magnitudes. However, the effective temperature and sur- 
face gravity derived from this model produced a poor fit 
to the UV spectrum. Then, the UV spectrum was fit 
independent of the photometry (i.e., T^s and log 5 were 
allowed to vary). This produced a much better fit with 
exception of two carbon lines at 1657 A and 1930 A (see 
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Fig. 7.— Spectral plot of WD 1009-184. The inset plot displays 
the spectrum (thin line) in the region to which the model {thick 
line) was fit (assuming a slight hydrogen abundance of log (H /He) 
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Figure 6). The asymmetry in the observed spectrum is 
not reproduced in the model fit likely because of the fail- 
ure of the impact approximation for the van der Waals 
broadening used in the model (Koester et al. 1982). With 
the Toff and \ogg fixed from the spectroscopic fit, the 
photometry was again fit to arrive at the final physical 
parameters found in Table 4, including an abundance 
of log (C/He) = -5.55 ± 0.12. We note that the dis- 
crepancy in the temperatures derived from the UV spec- 
trum and the optical/near- IR photometry is fairly small 
(-10,250 K vs. -11,300 K), though significant. It is pos- 
sible that the models fail to address some component of 
the input physics, such as a missing opacity source that 
gives rise to this discrepancy between UV and optical 
effective temperature determinations. 

WD 0821-669: a cool DA WD that was uncovered 
during a trawl of the SuperCOSMOS Sky Survey (SSS) 
database (Subasavage et al. 2005a, 2007). It is one of 
the oldest and nearest WDs of the new 25 pc members - 
6.22 ± 0.16 Gyr at a distance of 10.65 ± 0.12 pc. Thus, 
it is now the 23rd nearest WD system. 

WD 1009-184: a DZ WD that is difhcult to model 
accurately because its effective temperature is near the 
point at which additional pressure effects become impor- 
tant and are not included in the DZ models. Also, trace 
amounts of hydrogen may exist in the atmosphere but at 
levels too low to detect spectroscopically, yet affect the 
profile of the calcium absorption as discussed in Dufour 
et al. (2007). Our best fit using both photometry and 
spectroscopy produces a Tes = 5940 ± 280 K, including 
trace amounts of hydrogen [log (H/He) = —3] and an 
abundance of log (Ca/He) = -10.37 ± 0.20 (see Figure 
7). However, it is likely that additional pressure effects 
have an impact, so the physical parameters should be 
considered preliminary estimates. 

WD 1036-204: a peculiar DQ WD that has carbon 
Swan absorption bands; however, the absorption is sig- 
nificantly deeper in this object because of the presence 
of a large magnetic field (Bues 1999). Again, no atmo- 
spheric modeling was possible but this object will serve 
as another important test case for model revisions in the 
future. 

WD 2008-600: a DC WD that was recently discov- 
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Fig. 8. — Light curve of WD 2040—392 (top panel) normalized 
by its mean instrumental magnitude and the Fourier (amplitude) 
spectrum (bottom panel) identifying the dominant pulsational pe- 
riod. 



ered to have a flux deficiency in the infrared because of 
collisions by molecular hydrogen with helium, similar to 
WD 0038-226. Subasavage et al. (2007) have modeled 
this object using a mixed H/He composition and a pre- 
liminary trigonometric parallax. The updated physical 
parameters listed in Table 4 replace the preliminary par- 
allax with the more accurate one presented here, and we 
derive a helium abundance of log (He/H) = 2.60. Also, 
this object has the largest tangential velocity (112 km 
sec~^) of the systems presented here. 

WD 2040-392: a DA WD that is hsted in the 
McCook-Sion White Dwarf Catalog (McCook & Sion 
1999)^ but had no follow-up observations until Subasav- 
age et al. (2007) obtained a photometric distance esti- 
mate of 23.1 ± 4.0 pc (in excellent agreement with the 
trigonometric parallax distance of 22.63 ± 0.51 pc pre- 
sented in this work). In addition, as discussed in § 3.1, 
this object is variable at the —2% level. The physical 
parameters listed in Table 4 (i.e., T^g = 10,830 ± 310 
K, log g = 8.03 ± 0.04, mass = 0.62 ± 0.02 Mq, and 
absolute F = 11.98 ± 0.06) are entirely consistent with 
new pulsating ZZ Ceti WDs discovered by Gianninas et 
al. (2006). 

For confirmation of ZZ Ceti-type pulsations, data were 
acquired at the CTIO 0.9m using the same central quar- 
ter of the CCD as used for parallax data. Observations 
were taken in white light to maximize the signal for the 
target and eleven reference stars with a temporal resolu- 
tion of ~1 minute. Relative aperture photometry was 
performed on the target and reference stars using an 
aperture diameter of 12". The data span one hour, dur- 
ing which three cycles of a regular pulsation are clearly 
evident (sec Figure 8), thus confirming the target is a ZZ 
Ceti pulsator. The Fourier (amplitude) spectrum identi- 
fies the dominant period of ~980 seconds with an ampli- 
tude of -3%. 

WD 2251-070: a cool DZ WD for which the at- 
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http://heasarc.nasa.gov/W3Browe/all/mcksion.html. 
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mospheric model appropriate for DZs fails to reproduce 
the observed spectrum (sec Dufour et al. 2007) likely be- 
cause of additional pressure effects not accounted for in 
the model. Also, this object's effective temperature is at 
the limit of the model grid. Physical parameters should 
be considered preliminary estimates. 

WD 2336-079: a DA WD that Gianninas et al. 
(2006) recently discovered is a pulsating ZZ Ceti WD. 
Teff and logg listed in that publication are uncertain. 
With the high quality parallax presented here, the pa- 
rameters are now better constrained and in reasonable 
agreement with Gianninas et al. (2006). Also, with a 
minuscule proper motion of 33.6 ± 1.0 mas, this object 
is the slowest moving WD known in the 25 pc WD sam- 
ple. In fact, it was first cataloged by Giclas et al. (1975) 
and labeled as a suspected WD simply based on its blue 
color. Indeed, this object has the largest effective tem- 
perature of all of the systems presented here. Thus, this 
object provides additional support for the possibility that 
a sample of nearby WDs with little to no proper motions 
may exist, especially given that the majority of WDs in 
a volume-limited sample are cooler such that their colors 
alone are unexceptional. 

WD 2351-335: a DA WD that has both a M4V pri- 
mary (LHS 4039) and a recently discovered active M8.5V 
companion (APMPM J2354-3316C, Scholz et al. 2004). 
The WD has a separation of 6'.'5 at position angle 182.2° 
from the primary and the M8.5V has a separation of 
102'.'8 at position angle 91.3° from the primary. Par- 
allax data were taken in / to optimize observations for 
the M4V primary. Thus, the secondary WD was poorly 
exposed in the images giving rise to the largest paral- 
lax error of those presented here. The weighted mean 
parallax listed in Table 4, which represents our two mea- 
surements, should be taken as the parallax of the system. 
The tertiary M8.5V was unknown at the time parallax 
observations began (mid-2003) so, while the object is vis- 
ible in our frames, exposure times were not sufficient to 
obtain usable astrometry. 

6. DISCUSSION 

We have increased the number of WD systems with 
reliable trigonometric parallaxes within 25 pc from 99 
to 119 (20%). All are in the southern hemisphere (see 
Figure 2), thus, we are shrinking the disparity between 
the number of nearby WDs in the two hemispheres illus- 
trated in Table 2. In fact, this effort has nearly doubled 
the number of 25 pc WDs in the —30° to —90° quadrant 
of the sky. While 13 of the new 25 pc members have been 
known WDs for many years, seven are recent WD dis- 
coveries (Subasavage et al. 2007, 2008). Trigonometric 
parallax determinations are underway for an additional 
~30 WD systems estimated to be within 25 pc from those 
publications as well as others (e.g., Holberg et al. 2008), 
and it is likely that we will continue to populate the 25 
pc WD sample with reliable members. 

A number of interesting objects have been uncovered 
through the collection of parallax data presented here 
coupled with previous efforts' data. For instance, the 
parallax of WD 0121—429 places a considerable con- 
straint on its mass of 0.41 Mq assuming a single WD. 
Because this value is likely too small to have formed 
via single star evolution, the parallax gives additional 
weight to the hypothesis adopted from spectral analy- 



ses that this object is an unresolved double degenerate. 
Also, a new ZZ Ceti pulsating WD was identified (WD 
2040—392). Finally, the high quality parallaxes for two 
cool WDs that display continuum absorption in the near- 
IR (WD 0038-226 and WD 2008-600) show that trace 
amounts of hydrogen in a helium- dominated atmosphere 
is the most likely scenario to explain this poorly under- 
stood phenomenon. These objects and others like it will 
serve as empirical checks and permit revisions to atmo- 
spheric models for the coolest WDs. 

In the course of collecting long time-series astrometric 
data for nearby WDs, we hope to identify new systems 
that have astrometric perturbations from unseen com- 
panions. In particular, identification of new double de- 
generate systems that are resolvable with high resolution 
astrometric instruments (i.e., Hubble Space Telescope^s 
Fine Guidance Sensors) will permit accurate dynamical 
mass determinations for two WDs at once. To date, only 
three WD systems (Sirius B, Procyon B, and 40 Eri B, 
Provencal et al. 2002) have dynamical masses known to 
better than 5%. In two of these cases, the largest uncer- 
tainties in the mass error budgets are the distance deter- 
minations even though all three systems are among the 
50 nearest systems to the Sun. Also, these three systems 
serve as the only reliable empirical verifications of the 
theoretical WD mass-radius relation commonly used in 
WD modeling (e.g., WD ages via cooling models). Thus, 
accurate trigonometric parallaxes of nearby WDs are es- 
sential to finding additional systems whose masses can 
be well constrained. 

Note added in manuscript. — Two recent trigonomet- 
ric parallax works were published during the review pro- 
cess of this manuscript that are relevant to the 25 pc 
WD sample. Gatewood & Coban (2009) present trigono- 
metric parallaxes for 21 systems including four new WD 
systems within 25 pc (WD 0423-^044, WD 0511-1-079, 
WD 1309+853 - the previous parallax had a parallax 
error larger than 10%, and WD 2047-1-372). Lepine et 
al. (2009) present trigonometric parallaxes for 18 sys- 
tems including two new WD systems within 25 pc (WD 
1814-1-134 and WD 2322-M37). Thus, the 25 pc WD 
sample contains 105 systems prior to this work and 125 
systems including this work. 
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TABLE 5 

ASTROMETRIC RESULTS. 



WD ir(rel) ^(corr) 7r(abs) P.A. V^^^ 

Name R.A. (J2000.0) Decl. (J2000.0) Filter ^Vsca ^frm Coverage Years W^^f (rrras) (mas) (mas) (mas yr^^) (dog) (km s^^) Notes 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) 
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Object's previous trigonometric parallax of 101.20 ± 10.40 mas in YPC failed to meet the fractional parllax error constraint of ] 

New ASPENS member - within 15 pc. 

Affected by cracked V filter discussed in § 3.2.3. 

Common proper motion companion to WD 0628-020. 

Hvppnrcos parallax for companion LTT 2976 of 51.52 ± 1,46 mas. 

Hipparcos parallax for companion LHS 2231 of 58.20 ± 1.67 mas. 

Due to reddening of the reference stars, an average correction to absolute parallax was adopted (see § 3.2). 



Not affected by the cracked V filter discussed i 



'Commoi 
■■yPC pa 
^YPC pa 
^Commo: 
"^YPC pa 
"yPC pa 
°YPC pa 



proper motion companion to WD 2351- 
allax of 155.00 ± 2.10 mas. 
allax of 112.40 ± 2.70 n 



§ 3.2.3 because observations began after filter switch. 
-335. 



of 141.10 ± 8.40 ] 
of 112.70 ± 9.70 ] 
of 218.30 ± 6.70 ] 
^YPC parallax of 123.70 ± 4.30 i 
'^YPC parallax of 127.40 ± 6.80 i 



. WD 0738-172. 



and a Hipparcos parallax of 217.01 ± 2.40 mas. 



